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We generalize modern ideas about the duality between Wilson loops and scattering amplitudes in 
TV = 4 SYM to large- iV (or quenched) QCD. We show that the area-law behavior of asymptotically 
large Wilson loops is dual to the Regge-Veneziano behavior of scattering amplitudes at high energies 
and fixed momentum transfer, when quark mass is small and/or the number of particles is large. 
We elaborate on this duality for string theory in flat space, identifying the asymptotes of the disk 
amplitude and the Wilson loop of large- iV QCD. 
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This Letter is inspired by a remarkable recent discovery 
of the duality between Wilson loops and scattering ampli- 
tudes in AT = 4 super Yang-Mills (SYM) theory (see Q 
for a review of this subject). SYM differs from QCD 
by the contents of matter fields (6 scalars and 4 spinors 
in the adjoint representation of the SU(7V) color group, 
thereby providing an extended Af = 4 supersymmetry) 
and has attracted a significant interest over the last three 
decades as a toy model for certain aspects of QCD, in 
particular, for the relation between QCD and strings. 
Adding the extra fields makes the dynamics of SYM much 
simpler than that of QCD and it enjoys the famous anti- 
de-Sitter-space/conformal-field-theory (AdS/CFT) cor- 
respondence which, in particular, relates SYM Wilson 
loops to an open superstring in AdS5 S 5 [2\. While 
essential ingredients of QCD — asymptotic freedom and 
quark confinement — are not present in Af = 4 SYM, it 
captures many features of QCD perturbation theory. 

The (finite part of the) 4-gluon on-shell scattering am- 
plitude in SYM has the form 

A(s,t) = A tICC e fW^(s/t) (1) 

(where s and t are usual Mandelstam's variables) as was 
conjectured 0] on the basis of three-loop calculations. 
To explain Eq. ([I]), the Wilson- loop/scattering- amplitude 
(WL/SA) duality was introduced [4| at large 't Hooft cou- 
plings A, which has been then advocated in SYM pertur- 
bation theory 5]. This duality states that the scatter- 
ing amplitude (divided by the kinematical factor A tTee ) 
equals the Wilson loop for a rectangle whose vertices Xi 
are related to the momenta pi of scattering gluons by 

Pi = K (xi - Xi-i) , (2) 

where K — l/2ira' is the string tension. 

The function /(A) also appears in the anomalous di- 
mensions of cusped Wilson loops and operators of twist 
two. It has been recently found as a solution to the equa- 
tion derived from spin chains. Its perturbative solu- 
tion reproduces the three 0, 0| and four || loop SYM 



results, while numerical [9( and analytical [lfj solutions 
reproduce the VA behavior of /(A) for large A originally 
found [ll| using the AdS/CFT correspondence. The next 
orders in 1/ \f\ also agree with the superstring calcu- 
lations [12] thereby providing a remarkable test of the 
AdS/CFT correspondence. 

Our goal in this Letter is to find out what features 
of the described WL/SA duality (if any) remain valid 
for QCD and, in particular, how is it possible to main- 
tain the relation of the type © which would relate large 
momenta in scattering amplitudes with loops of large 
size. Of course this is not possible in QCD perturba- 
tion theory, where \p\ ~ l/\x\ because of dimensional 
ground. But nonperturbatively a dimensional parameter 
K « (400 MeV) 2 appears in QCD, which shows up in the 
area-law behavior of asymptotically large Wilson loops: 

W(C) UIS £ C e ~ KS ^ c \ (3) 

where <5 , m in(C) is the area of the minimal surface bounded 
by C, that results in confinement. Strictly speaking, this 
requires large N or the quenched approximation. 

As is well-known by now, a string theory, which QCD 
is supposedly equivalent to, is not the simplest Nambu- 
Goto string. Some extra degrees of freedom living on the 
string are required which are most probably conveniently 
described by a presence of extra dimensions. The asymp- 
totic behavior ((3]) is nevertheless universal for large loops. 
Also, there is a considerable amount of evidence from 
lattice gauge calculations in 2+1 and 3+1 dimensions 
for various N that the Nambu-Goto action describes the 
behavior of the Wilson loops quite well and the transi- 
tion from perturbative to stringy behavior takes place 
"at surprisingly small distances" [13j |. There also exists 
a number of other comparisons between results from the 
Nambu-Goto action, e.g. between the closed string spec- 
trum, and SU(A^) for various N (see [3] and references 
therein). This action has the well-known anomaly for 
d 7^ 26, which however is suppressed for long strings [HI]. 
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This remarkable success of the Nambu-Goto string in 
flat space as an effective action leads us to reconsider the 
relation between the Wilson loop W(C) and the corre- 
sponding string wave functional. We then obtain scatter- 
ing amplitudes in large- N (or quenched) QCD by prop- 
erly summing W(C) over paths and find that the WL/SA 
duality holds in a kinematical region of large s and fixed 
t when only large loops, for which the area law ([3]) sets 
in, are essential in the sum over paths. Thus obtained 
scattering amplitudes, quite involved in general, are of 
the Regge-Veneziano type when the quark mass is small 
and/or the number of external particles is large. 

Our starting point is the standard representation of 
Green's functions of M colorless composite quark opera- 
tors (e.g. q{xi)q{xi)) in terms of the sum over all Wilson 
loops passing via the points Xi (i = 1, . . . , M), where the 
operators are inserted: 



M 



g= m^M a 

\i=l 

„r Af-a 



dT e- mT 



[ dr M -i TT / *dr, f Vz{t) J[z(t)] W[z(t)}. (4) 
Jo Jo J 



«(O)=*(T)=s> 
»(T 4 )= s . i 



Here the weight for the path integration is 



J[z{t)} = / Vk{r) sp Pe 



ij^dr [i(r)-fc(r)- 7 (T)-fe(T)] 



(5) 



for spinor quarks and scalar operators. In Eq. (j4|) VF(C) 
is the Wilson loop in pure Yang-Mills theory at large N 
(or quenched), m is the quark mass and r is the proper- 
time variable. For finite N, correlators of several Wilson 
loops have to be taken into account. The derivation of 
this formula and the references can be found in [l6j |. 

The on-shell M-particle scattering amplitudes can be 
obtained from the Green function (01 by the standard 
Lehman-Symanzik-Zimmerman reduction. When mak- 
ing the Fourier transformation, it is convenient to repre- 
sent M momenta of the (all incoming) particles by the 
differences Api = pi-\ — pi . Then momentum conserva- 
tion is automatic while an (infinite) volume V is pro- 
duced, say, by integration over xq. It is convenient to 
introduce a momentum-space loop P^(t) which is piece- 
wise constant: 



p(r) = pi for T t < T < Tm 



(6) 



Because p(r) = - £\ ApiS(T - n) with Api = p^i - p 4 , 
we write in the Fourier transformation: Api ■ Xi = 
J drp(r) • z{t) which is manifestly parametric invariant. 
Making the Fourier transformation, we obtain 

G (Api , . . . , ApM ) = / dTe- mT [ dr M -! 

Jo Jo 

xf[ f^&n f Vz(j) e 1 ^ 'dri(r).p(r) J[ Z ( T )]W[z{t)], 
z(0)=z(T)=0 (7) 



where p(r) is piecewise constant as in Eq. ([6]) . We do not 
integrate over z(0) = z(T) which would produce the (in- 
finite) volume factor because of translational invariance. 

To calculate the scattering amplitudes, we have to sub- 
stitute the area-law behavior © of asymptotically large 
Wilson loops into Eq. ([7]) and to integrate over the paths. 
In general, this would lead us to very complicated inte- 
grals but the calculation drastically simplifies if to use 
the representation of the minimal area as a boundary 
functional that was introduced by Douglas in his 
celebrated solution of the Plateau problem. We shall use 
one of the equivalent forms of the Douglas functional: 



I f T 

A \°\ = ~T~ / dT 1 dT 2 x(T 1 ) ■ x(t 2 ) 
^ Jo 



x In (1 - cos [2tt Hn) - a(r 2 )) /T]) , (8) 

where < <t(t) < T is a reparametrization (cr'(r) > 0). 
The functional is to be minimized with respect to 
er(r) with the minimizing function <t*(t) being of course 
contour-dependent. Then A[ov] is equal to the minimal 
area S m m(C), while in general A [a] > A[cr*\ = 5 m in(C). 

In fact (jlj is well-known as the classical boundary 
action in string theory. It appears for the tree-level 
disk amplitude with Dirichlet boundary conditions in the 
Polyakov string formulation after integrating over the 
string fluctuations inside the disk, i.e. over X(r, 9) with 
r < 1, < 9 < 2tt, and fixing the value X(l,9) = x{9) 
at the boundary. The appearance of the function a* (9) 
is related to a subtlety associated with fixing conformal 
gauge [IH . The decoupling of the Liouville field is pos- 
sible only in the interior of the disk, while its bound- 
ary value determines the function cr*(0) at the classi- 
cal level. The path integral over the boundary value 
of the Liouville field then restores the invariance under 
reparametrizations of the boundary in quantum theory. 

Motivated by this fact, Polyakov [l9j proposed to iden- 
tify the Wilson loop in large-iV QCD with the tree-level 
string disk amplitude integrated over reparametrizations 
of the boundary contour. It is convenient to confor- 
mally map the disk into the upper half-plane, so the disk 
boundary is mapped into the real axis parametrized by 
t(r) = tan(7rr/T), ~oo < t < +oo. Then we write 

W{C) = J Vs{t) exp ( £ J + dtidt a ±(ti)-x(ta) 

xlnKtO-sfc)!), (9) 

where the path integral over s(t) (with s'(t) > 0) restores 
the invariance under reparametrizations. 

In spite of the fact that the right-hand side of Eq. ([9]) 
is derivable for bosonic string in d = 26 or superstring in 
d = 1 0, we shall use it only for asymptotically large loops 
or, equivalently, very large K, when the integral over 
reparametrizations has a saddle point at s(t) = s*(t). 
This is crucial for reproducing Eq. ([5]). 
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It is easy to calculate a (reparametrization-invariant) 
functional Fourier transformation 



W\p(-)] = / Vx e ^ p - dx W[x(-)] 



(10) 



of the disk amplitude © for piecewise constant p(t). 
Substituting ^ into Eq. (fTU|) and performing the Gaus- 
sian integration, we get 



W[p(-)} 



+ 00 



Vs(t)exp(a' dt^pih) -p(t 2 ) 

J — OO 

xln|«(*i)-«(*2)f) (11) 



which is of the same form as ((9]) only with K replaced by 
l/K = 2W. 

Since p(t) — pj at the j-th interval for the stepwise 
discretization, the only effect of the reparametrization is 
to change the values of tj 's for Sj 's keeping their cyclic 
order. This is a discrete version of the reparametrization 
transformation. Note that the stepwise discretization of 
x(t) itself is not possible since it would violate the conti- 
nuity of the world-line of the string end. 

The stepwise discretization ([6]) naturally results in the 
M -particle (off-shell) Koba-Nielsen amplitudes which are 
invariant under the SL(2;R) projective transformation 
s (as + b)/ (cs + d) with ad — be — 1 because the pro- 
jective group is a subgroup of reparametrization trans- 
formations. To derive them, we first note that 



dtidt 2 p(*i)-p(*2)ln|»(*i)-*(*a)| 

1 r +o ° d Sl ds 2 



[p(t( Sl ))- P (t( S2 ))r (12) 



for the integral in the exponent in (|11[) . The integration 
over Si or s 2 on the right-hand side has divergences when 
they lie on adjacent sides k = I ± 1. If we omit the sides 
with k — I ± 1, then the integrations over s\ and s 2 are 
perfectly finite resulting in 



1 



Sk A r A (Pk-Pl) 
dsi / ds 2 



(si- S2 ) 2 



= ^2 A Pk ■ Api log \sk - si\ 

k^l±l 

( Sj - Sj-i)(s j+ i - Sj) 



+E A ^ 2l °g- 



(s j+1 - s 3 _i) 



(13) 



which is projective invariant. 
Choosing the measure to be 



Vs 



[[d Sl 



Si+i - Si-i) 



< (14) 
(Si - Si-i)(s i+ i - Si) 

which is also invariant under the projective transforma- 
tion, we arrive at 



W(A Pl ,...,A PM )= f X\d Sl \[\ Sk - Sl \ 



a'Ap k Api 



n 



(Sj - Sj-iXgj+i - Sj) 

(Sj+i - Sj-!) 



a'Ap 2 ,-! 



, (15) 



where the integration over Sj emerges from the path 
integral over reparametrizations in Eq. (Ill[) . This is 
known as the Lovelace choice j2(| (see |2l|). that re- 
produces some projective-invariant off-shell string ampli- 
tudes known since late 1960's. The more familiar on-shell 
tachyon amplitudes can be obtained from Eq. (TT5]) by set- 
ting a'Ap 2 = 1. 

Fixing in Eq. (|15p the remaining SL(2;K) invariancc 
in the standard way, we obtain the scalar amplitudes in 
the Koba-Nielsen variables. For the case of 4 scalars this 
reproduces the Veneziano amplitude 

A(A Pl , Ap 2 , A P3 , A P4 ) = / dxx- a ^-\l - x)- a ®-\ 



where a(s) = a's + 1 and 



(16) 

= -(A Pl + Ap 2 ) 2 , 
t = — (Ap 2 + Ap 3 ) 2 are usual Mandelstam's variables 
(for Euclidean metric). Here the tachyonic condition 
a'Apj = 1 has not to be imposed. While Eq. (TT5|) re- 
sults in a(0) = 1, an arbitrary value of the intercept a(0) 
can be reached by properly changing the measure (|14| . 

We are now in a position to perform the main task 
of this Letter: to substitute the area-law behavior © 
of W(C) into the path integral ([7]) and to find out for 
what momenta the asymptotically large loops dominate. 
As we shall see, typical momenta will be large for large 
loops. As is already explained, we substitute Eq. ((3]) by 
Eq. © which gives the same for large loops (or large 
K = 1/2W). 

Interchanging the order of integration over z(t) and 
<j(t) (or s(t)) and easily doing a Gaussian path integral, 
we obtain 



G (A Pl , . . . , Ap u ) = dTe- mT dr M -i [] jf *" dr 4 J Vo(r) J Vk(r) 

X Sp P e^^ Tdr i/o Tdr 2(' ;: ( T i)+P( T i))-(fe( T 2)+p(^))ln(l-cos[^( ( T(T 1 )- ( T(r 2 ))])-i/ T dr7(r)-fc(T)^ Qys 



4 



This expression is rather close to the disk amplitude (fTTj) . 
except for the additional integration over k. But for the 
case where m is small (or M is very large), the integral 
over T in Eq. (fTT|) is dominated by T ~ (M — l)/m which 
is large for m small. This is because JX 1 dr^ ~ T 711-1 . 
Noting that typical values of k ~ 1/T are essential in the 
path integral over k for large T, we can disregard fc(r) in 
the exponent in Eq. (jTTJ) so the integral over k factorizes. 
Making the change of the variables from a to s, we finally 
obtain from Eq. (|17|) the product of the momentum-space 
disk amplitude times factors which do not depend on 
p. Therefore, Eq. (|17p exactly reproduces for piecewise 
constant p(r) the (off-shell) Koba-Nielsen amplitude (fT5|) 
as m — > 0! 

It still remains to discuss in what kinematical region 
of momenta Api our derivation is legible, that is only 
asymptotically large loops are essential in the path inte- 
gral over z(t) in Eq. 0. A physical intuition suggests, 
from the spectrum of a classical string, this should be the 
case at least for asymptotically large s and large t < s. 
This indeed agrees with our formulas, where the value 
of ev(0) is not essential in this region. But the domain 
of applicability of our approach is broader and extends 
to large negative values of t. However, when — t <C s be- 
comes large enough, there are no longer reasons to expect 
the contribution of large loops to dominate over pertur- 
bation theory, which comes from integration over small 
loops in Eq. (|7|). Therefore, our formula for the 4-point 
scattering amplitude is valid only for asymptotically large 
s and fixed t {\t\ > 1/a 1 ), associated with small angle or 
fixed momentum transfer. The tachyon issue, which is a 
short distance phenomenon [ljj . is then irrelevant. For 
smaller values of |t| < 1/a' the results become sensitive 
to the choice of the measure in the ansatz ©. 

As distinct from previous approaches to reggeization 
in perturbative QCD, in particular from that based on 
the evolution equation [22j] for Regge trajectories, our 
approach deals with large loops usually associated with 
nonperturbative effects. Actually we are dealing with the 
quark-antiquark Reg ge t rajectory, whose QCD calcula- 
tion was pioneered in [23( , rather than with the Pomeron. 

When m is not small and/or M is not large, one should 
consider the full expression (fT7|) . We can split there 
the k— integral into two domains with small and large 
k. Then the former will appear as a Regge- Veneziano 
behaved factor coupled to the rest of the integrand. 

Thus, in conclusion we see that the area-law behavior 
of Wilson loops is dual to the Regge-Veneziano behav- 
ior of scattering amplitudes at high energies and small 
angles, when quark mass is small and/or the number of 
produced particles is large, but this ceases to be valid 
when the momentum transfer is large. This is how the 
exponential falloff of the 4-particle amplitude with large 
—t ~ s, which is unavoidable in string theory [24|, does 
not happen in our consideration. 
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